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Abstract Performance degradation remains as one of the primary limitations for 
practical applications of proton exchange membrane (PEM) fuel cells. The perfor- 
mance of a PEM fuel cell stack is affected by many internal and external factors, such 
as fuel cell design and assembly, degradation of materials, operational conditions, and 
impurities or contaminants. Performance degradation is unavoidable, but the degrada- 
tion rate can be minimized through a comprehensive understanding of degradation and 
failure mechanisms. In present work, a single PEM fuel cell for stationary applications 
is investigated. Membrane and catalyst layers (anode and cathode electrodes) are con- 
sidered as critical components that affect the degradation of the cell. The model used in 
this work diagnoses degradation of MEA (platinum degradation in catalyst layers and 
membrane thinning and dehydration in polymer membrane), and by considering the 
degradation over operating time, estimates power density over system lifetime. In this 
paper, also three optimization model with different objective functions are developed 
tomaximize total energy production. The results show that by continuously optimiz- 
ing the operating conditions, total energy generation of the system will increase up to 
3.16 96. 
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1 Introduction 


Extending fuel cell total energy (kw.hr energy produced during fuel cell lifetime) is 
a necessary objective for reducing fuel cell per unit electricity cost. Total generated 
energy is increased by extending life time and operating at maximum possible output 
power. Fuel cell life time and output power are affected by degradation mechanisms. 
Reducing the rate of degradation mechanisms can be achieved by implementing an 
appropriate optimization strategy that prevents the system from being placed in the 
operating regions that cause failure. To obtain the optimization strategy, first a degra- 
dation model is developed. The degradation procedure of fuel cell stack components 
and its impact on electricity generation through the remaining lifetime, in most cases, 
has been insufficiently understood. The lengthy required testing time, lack of under- 
standing of most degradation mechanisms, the difficulty of performing in-situ, and 
non-destructive structural evaluation of key components make the issue difficult to 
investigate. The fuel cell membrane electrode-assembly (MEA) plays a vital role in 
the overall lifetime achieved by a stack [1]. MEA is the main component of our study. 
Our approach to improve PEM fuel cell total energy production is to investigate the 
degradation mechanisms and optimizing operating conditions which maximize cell 
total energy regarding degradation mechanisms. The MEA performance shows degra- 
dation over operating time, which is dependent upon fuel cell design and assembly, 
degradation of materials, operational conditions and impurities or contaminants. To 
model the degradation mechanisms, itis necessary to understand degradation processes 
of the PEMFC components. Various studies have been conducted in this field. Their 
differences are mostly in studied components, experimental or model based degra- 
dation diagnosis and the purpose of degradation diagnosis. Borup et al. [2] made a 
review on PEM fuel cell components (membrane, electrodes and bipolar plates) degra- 
dation. In this review, some of the targets for PEM fuel cell systems, as indicated by 
the US DOE (Department of Energy) and Japanese NEDO (New Energy & Indus- 
trial Technology Development Organization) are presented showing the importance 
of PEM fuel cell degradation research. The purpose of their paper is to investigate the 
effects of operating conditions (Impurity and Subfreezing) on fuel cell durability using 
experimental data of durability tests. Yousfi et al. [3] reviewed PEM fuel cell catalyst 
degradation using available experimental research. In their paper, the effect of fuel 
and oxidant starvation on catalyst and carbon support degradation is investigated and 
some characterization methods or mitigation strategies for catalyst layer are presented 
to build catalyst illustrative fault tree. Darling et al. [4] developed a mathematical 
model of oxidation and dissolution of supported platinum catalysts in PEM fuel cells 
to study electrode potential and particle size effects on catalyst stability. Chen et al. 
[5] studied membrane degradation. In their study, experimental degradation data are 
collected, the membrane samples are analyzed, and the effects of Fe^* concentration 
and temperature on membrane degradation are discussed. In their dissertation, both 
experimental and modeling methods are applied to investigate the chemical durability 
of PFSA membranes for fuel cell applications. Bruijn et al. [6] reviewed the degrada- 
tion mechanisms of membranes, electrodes, bipolar plates and seals. In their study, by 
collecting long-term experiments, the relation between degradation mechanisms and 
the operating conditions (voltage, pressure, humidity and temperature) was discussed. 
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Schmittinger et al. [7] presented the effect of humidity on fuel cell components (catalyst 
layers, the gas diffusion media and the membrane) degradation in long term perfor- 
mance. The proposed research shows that the main cause of short life and performance 
degradation is poor water management. Jayson et al. [8] optimized peroxide mitigation 
strategies, mechanically reinforced membrane and chemically stabilized membrane 
by decreasing cost and performance impact. They try to improve the lifetime of PEM 
for fuel cell by seeking technologies that will prevent membrane failure. 

In the present work, unlike the previous studies, membrane and catalyst degradation 
that are dominant, especially in stationary applications are modeled simultaneously, 
and the effects of temperature, humidity, and pressure on degradation mechanisms are 
investigated. 

In most researches, optimum technologies and structure parameters regarding cost 
and performance degradation are obtained whilein this paper structure parameters 
has been fixed and the effect of operational conditions on degradation mechanisms is 
discussed. To do this, a model has been developed which includes membrane chemical 
degradation, platinum growth in catalyst layer and electrode reduction in active surface 
area. Using degradation model, voltage and current are estimated as a function of 
catalyst and membrane degradation. Degradation itself is expressed as a function of 
operating conditions. Afterwards, by optimizing with three objective functions the 
optimal operation conditions which maximize power output through fuel cell life time 
are presented. 

In short, the main innovations in this paper include: 


e A dynamic model which considers the decreases in catalyst active surface area 
and the membrane thinning simultaneously 

e The effect of membrane and catalyst degradation on fuel cell performance 

e Optimized operational conditions to maximize the total energy over fuel cell life 
time regarding degradation mechanisms 

e History based optimization strategy that maximize system total energy. 


2 Mathematical modeling 
2.] Performance modeling 
The general performance model which is used in most researches to estimate fuel cell 
performance is presented in Eqs. 1—2 [9]. In this model degradation are not consid- 
ered. Therefore, equations are time independent. The fuel cell output power density 
is calculated as Eq. 1. 

P = Vfce xi (1) 
where output voltage (V ) of PEM fuel cell is greatly affected by Ohmic, concentration, 
and activation polarizations. The operating cell potential can be formally expressed as 


Eq. 2. Related terms are presented in Table 1 [10]. 


Vee = Vrev — Vohm — Vact — Veone (2) 
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Table 1 Voltage term equations 


Voltage terms Equation 











Ves 1.229 + (4.308 x 1075)7 x In | Pry (Poy) | — (8.453 x 10-4)(T — 298.15) 3) 
181.6] 1+0.03i-+0.062( 45,225] 
(.—0.634—3i) exp] 4.18( 393) | 
RT VQ iu. 
Veross&act onF In ig (5) 
Vconc RE (1+ 3) #5 (6) 


2.2 Degraded performance modeling 


As fuel cell works, due to chemical reactions, degradation occurs. Degradation affects 
fuel cell performance. In this section, anew model is developed which considers degra- 
dation effects on fuel cell performance. Considering degradation makes equations time 
dependent. 

Recognition of degradation mechanisms and their dependence on operation condi- 
tion is critical for understanding fuel cell longterm performance. In order to determine 
power density through time, the value of the voltage and current regarding degradation 
should be estimated as the function of time. Power density through fuel cell life time 
can be expressed as Eq. 7. 


P(t) = Vec(t) x i(t) (7) 
2.2.1 Voltage degradation 


In degraded performance model, overpotentials are time dependent and can be 
expressed in following equations. 


Vrc(t) = Vrey — Vonm(t) — Veross&act (t) — Veone — VECSA(t) (8) 


Degradation model of over potentials is developed based on components degradation; 
all losses are functions of the catalyst and membrane degradation. 

Ohmic loss Ohmic loss in PEM fuel cell is caused by resistance to hydrogen ions 
transition and 1s calculated as follow [11 |: 


Vonm(t) = i X rion(t) (9) 


181.6 [1 + 0.03i + 0.062(57;)2i?? ] 


| : (10) 
(à — 0.634 — 3i) exp [4.18(7—299) | 


l'ion(t) = tn(t) 


Ionic resistance is totally dependent on membrane water content and membrane 
thickness. In Eq. 10, A, represent membrane water contentwhich shows number of 
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water molecules with respect to SO?- molecules in membrane [12]. Acan be related 
to membrane relative humidity (RH) using Eq. 11[13]. 


à = (0.043 + 17.81 RH — 39.85R H? + 36.0R H?) (11) 


Here, is assumed that membrane relative humidity is equal to average of anode and 
cathode relative humidity. Relative humidity is estimated using Eq. 12. 

P, 

P sat 


RH = 





(12) 


Polymer electrolyte membrane thickness decreases during fuel cell life time. Different 
phenomena can cause membrane thinning. The most common reason is chemical 
degradation. It has been proposed that carboxylic end groups left over from the Nafion 
manufacturing process may be susceptible to attack by radical species generated during 
fuel cell reactions. The proposed mechanisms are as follows: 


Step 1: R-CF? COOH + OH: — R-CF»- + CO» + H2O 
step 2: R-CF, + OH. — R-CF20H — R-COF + HF 
Step 3: R-COF + H20 — R-COOH + HF 


Hydrogen peroxide which is an intermediate of the electrochemical oxygen reduc- 
tion reaction forms the radical species, such as hydroxyl radicals, also hydrogen or 
oxygen permeating to the anode catalyst layer may react and produce peroxide species. 
Produced radicals cause membrane thinning, release of fluoride ions, increased gas 
crossoverand consequently voltage degradation. 

Membrane thickness as the function of time through fuel cell lifetime is presented 
in Eq. 13 [14]. 


NC (5m + i) f(t) (13) 


Fraction of remaining electrolyte (f(t)) is proportional to hydrogen cross over 
according to Eq. 14 [14]. 


_ df(t) 
dt 





= Ki x Nm (t) x f(t) (14) 


Hydrogen crosses over current is related to membrane permeability and hydrogen- 
partial pressure as Eq. 16 [15]. 





in, (D) = 2F Np, (t) (15) 
“Cross = Py, Ap a 


Finally, hydrogen partial pressure and membrane permeability is determined 
according to Eqs. 17 and 18 [15]. 
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pe = Pci Po Dal d 1 
fy = (r) (ee (17) 
i Ig 2 

Pm m = 1.39 x 10 "e ??80/T (18) 


According to above explanation, the dependency of Ohmic loss to temperature, 
relative humidity, cell current and fraction of remaining electrolyte is expressed by 
Eq. 19. 


181.6i x (0.559 + 1) f (1) [1 + 0.03i + 0.062(45:)2i?? ] 
((0.043 + 17.81R H — 39.85 R H? + 36.0R H?) — 0.634 — 3i) exp [4.18( 15292) ] 
(19) 


Vohm (t) = 


Crossover and activation loss As any polymer membrane, the polymer electrolyte 
of a fuel cell allows some gas to permeate across it. Crossover rates of different gases 
have been measured for Nafion membranes and can be broadly described by Fick’s 
law. Hydrogen crossover to cathode incurs accelerating degradation mechanisms; in 
addition its effect on performance loss is added to activation losses, as shown in Eq. 
20 [16,17]. 


RT i + per oss (1) 
Veross&act (t) = —— In (an) (20) 
an F 10 


The correlation between Veross&gact(t)and operations conditions is presented as 
Eg: 21. 


RP. di do .f 139x076? NPA BL P 
Veross&act (t) = —— |n — Xft 1,0 a 2 
anF Vio (519 + D f (t) Pise: 


(21) 





Electrochemical active surface area loss One critical issue that affects voltage loss 
during PEM fuel cells operating lifetime is the electrochemical active surface area 
(ECSA) loss of carbon supported platinum nanoparticles at the cathode. The decrease 
in the electrochemical active surface area (ECS A) is corresponded to voltage loss by 
the Eq. 22. ao and a (t) are the ECSA at the beginning and at time t of fuel cell lifetime, 
respectively [6, 18]. 


V RT ao 
ECSA(t) = onF In p (22) 


Considering constant platinum (Pt) mass in electrodes, the correlation between Pt 
diameter and electrochemical active surface area is expressed by Eq. 23. ro and r (t) 
represent the particle initial radius and the particle radius at time t respectively [5]: 


%AECSA = m) — Sm (23) 
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Platinum radius changes mainly because of Pt dissolution and deposition in the 
catalysts. Four reactions lead to Pt dissolution and deposition. In the present study, 
active surface area degradation is modeled considering following reactions [19]: 


e Reaction 1: Pt electrochemical dissolution 
Pt — Pt?* + 2e- (U? = 1.188 V) 
e Reaction 2: Pt oxidation 
Pt + H20 > PtO + 2H* + 2e7 (US = 0.98 V) 


e Reaction 3: Diffusion of Pt ions and subsequent precipitation by hydrogen reduc- 
tion 


Pt +H —> Pie 
e Reaction 4: Chemical dissolution of Pt oxide 
PtO + 2H* — P+ + HO 


Last two reactions are not considered due to small rate of reaction 3 and the unclear 
kinetics of reaction 4. What is more, inclusion of proposed reactions does not affect 
the results measurably due to their slow kinetics [2]. 

According to Pt dissolution and deposition reactions rate, Pt particle radius pre- 
sented as Eq. 24 [2]: 


dr (t) M 
" —— (rı (t) + ra(t)) (24) 
t p 





In Eq. 24, r1(f) and r2(t) represent the rate of first and second reactions. rj (t)is 
calculated as Eq. 25 [15]. 


Qa 1 n1 F 
RT 


C pt2+ oc 1 n1 F 
(mes (-“ g == 22] 


where 0p;o(t) 1s the fraction of platinum surface covered by PtO and is obtained as 
Eq. 26. 


ri(t) = ky(1 — 0pio(t)) [exp ( ($1 — $2 — uso) (25) 


(26) 








dOpro(t) _ (22) B (220) dr (t) 
dt r(t) dt 


L max 


U(t) is the standard equilibrium potential of reaction 1 [2]. 
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o Abpr(t) 
pugne e 


LA(T) m U. 
1 (4) AF 


(27) 


Standard equilibrium potential is shifted away from the literature value for bulk plat- 
inum by the factor Ajtp; (t). The proposed factor accounts for the effect of surface 
tension of the platinum crystallite on standard equilibrium potential. 


Op Mp 
Anni) = pp, Q8) 
t 


The rate of second reaction is formulated as Eqs. 29-31, [2]. 


0 » 
r(t) = ka [exp (-220) exp (aa Hs Gio Use) 


2 
-onoo ( z “Ht o» C5 d = a vz) | (29) 
Cua ref 


Aupio() App;z(t) 





UD = UE US dd 
opto Mpio 
Anpil) = Alps t nono (31) 
t 


Active surface area loss during fuel cell life time is a function of temperature and 
Pt particles radius through time as is expressed in Eq. 32. Pt particles radius depends 
on its initial radius, cell potential and electrochemical reactions rate. 


= RT r(t) 
Vecsa(t) = —— In{ — (32) 


œan F ro 


Therefore, the degraded output voltage can be calculated by Eqs. (3, 6, 15, 21, 32). 
2.2.2 Current kinetic model 


During fuel cell operation, active surface area and hydrogen cross over current density 
increase. Therefore, net accessible current density during time, decreases as Eq. 33. 
In other words, if a fixed output current is desired, during operation, supply mass flow 
rate should increase [20]. 


i(t) =i x (1 — @pro(t)) — ig? (t) (33) 
Using equations in section (a), current density through time can be estimated as a 


function of temperature, pressure, fraction of remaining electrolyte and fraction of 
platinum surface covered by PtO (Eq. 34). 
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= —2980 
io cix a Broth) - (E EE (Pcr) (Eie t2) 





(510 4- 1) f (t) p 2 


inlet 


(34) 


3 Optimization 


The optimization objective is to maximize total energy production regarding degrada- 
tion mechanisms. Optimization model find the operating conditions which minimize 
the per unit energy cost. The decision variables which affect the power generation 
are cell temperature, anode and cathode relative humidity, as well as anode and cath- 
ode fuel pressure. The optimization model is also subject to several constraints. It 
is assumed that PEM fuel cell is designed to work in the specific temperature and 
pressure ranges. The proposed constraints are expressed in Eqs. 35-39. 


Tmin € T < Tmax (35) 

Panode,Min € Panode < Panode,Max (36) 
Pcathode,Min € Peathode € Pcathode,Max (37) 

R Acathode,Min € RHeathode < R Acathode,Max (38) 
R Hanode,Min € R Hanode € R Hanode,Max (39) 


Three objective functions are presented to compare the effect of controlling opera- 
tional conditions on fuel cell power generation and degradation mechanisms. In each 
proposed strategy, decision variables and constraints are the same but optimization 
objective functions are different. In this study, the objective function is based on both 
degradation and electrochemical phenomena; therefore, the effect of increasing of 
performance and cell degradation is estimated simultaneously. In the other words, no 
weighted methods are needed in our optimization formulation. 


3.1 Objective function 1 
Here, power density is maximized at the initial time (tọ) of fuel cell operation. In other 
words, degradation mechanism rate and power losses due to the degradation have been 


ignored. The operational conditions are constant through fuel cell lifetime. 


Max P (tg) = Max V (to) x I (to) (40) 
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3.2 Objective function 2 


To determine the operational conditions that maximize system total energy in 10,000 h 
of fuel cell lifetime, the second objective function is applied on the developed model as 
Eq. 41. In this strategy, the optimum operational conditions are obtained at the initial 
time and they are constant through fuel cell lifetime. 


110,000 110,000 
Max / P(t) -dt = Max J V(t) x I(t)-dt (41) 
to 10 


3.3 Objective function 3 


In this strategy, power density is maximized in every time step of fuel cell lifetime. 
It is clear that the operational conditions may be changed continuously through time. 
Equation 42 presents the third objective function: 


Max P(t) = MaxV (t) x I(t) (42) 


3.4 Optimization procedure 


Figure | points out the procedure of system optimization regarding degradation. The 
degradation rates are considered to determine the optimum values of pressure, tem- 
perature and humidity through fuel cell life time. Genetic algorithm is used and the 
optimum operation conditions for each objective function are obtained. 


4 Validation 


To solve a system of Eqs. 1—34, the fourth order Runge—Kutta method is used and the 
system variables are found. For validation, it should be noted that at present, the voltage 
degradation curve is the only experimental piece of data available in the literature that 
considers implicit degradation mechanisms. Therefore, voltage degradation curve is 
commonly used to validate the agreement between the model and experiments. Model 
results are verified with experimental results (at the operating conditions T = 348 K, 
Panode = Peathode = 1 bar, RHanode = RH cathode = 70 96) provided by S. Higashiguchi et 
al. [21]. According to Fig. 2, degradation curve of cell potential versus long runtime 
for the model is in good agreement with the experimental curve. For validation output 
power through time regarding degradation, Fowler [22]presented power density at 
discrete points of fuel cell lifetime (at the operating conditions T = 353 K, Panode = 
Pcathode = 3bar, RHanode = RH cathode = 100 %), the comparison is also presented in 
Table 2. As can be seen, the model results are nearly close to the experimental data 
provided by Fowler [22]. 
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— Model output 






0 2000 4000 6000 8000 10000 
Hour 
Fig. 2 Comparison between the modeling and experimental results for output voltage in PEM fuel cell life 


time 


Table 2 Comparison between the modeling and experimental results for power density degradation in 
PEM fuel cell 


Time (h) Model output Power density decrement (46) Power density decrement (96) 
(power density W/cm?) (model output) (experimental date) 

0 0.765 0 0 

1,500 0.6837 11 14 

3,000 0.6454 16 19 

4,500 0.612 2l 29 

6,000 0.5833 23 28 


5 Results and discussion 


To investigate the effects of operation conditions on cell performance loss due to 
system degradation, the proposed model is applied to the case study (T = 348 K, 
Panode = Peathode = 1 bar, RHanode = RH cathode = 70 96). Figures 3 and 4 present the 
main output of degradation model for the case study. Figure 3 presents the increment of 
degradation mechanisms. Three important observations stand out from Fig. 3: first, at 
the beginning of fuel cell operation lifetime the loss of active surface area is dominant; 
second, membrane degradation rate increases linearly during fuel cell lifetime; third, 
the catalyst degradation is almost more significant through fuel cell life time. 

Figure 4 indicates the behavior of output voltage and power density through the 
fuel cell lifetime. The rates of voltage and power density degradation are almost the 
same through fuel cell life time. It is also found that, after 10,000 h of operation, 22 96 
power density and approximately 9 96 output voltage are decreased. 

The base case parameters are presented in Table 3. Figure 5 presents the optimum 
temperature for proposed objective functions through the time. For the first objec- 
tive function which considers the maximum power density regardless degradation 
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Fig. 3 Degradation rate in PEM fuel cell lifetime 


100 voltage degradation (96) 
= =e Power density degradation (96) 
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=) 
= 


80 | 





75 | 
0 2000 4000 6000 8000 10000 


Time (hr) 


Fig.4 Voltage and power losses in PEM fuel cell lifetime 


mechanisms the maximum operating temperature is resulted. However, for second 
objective function the optimum temperature is decreased to minimize degradation 
rate which results in increase in total energy output of fuel cell. For third objective 
function the optimum temperature is determined for any desired time through fuel cell 
life time. The results show that the optimum operating temperature is decreased as 
fuel cell suffers degradation through time. 
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Table 3 Value of constants in 


optimization model ee a a 
TMax 373 K 
TMin 313 K 
Panode,Max 3 bar 
Panode,Min 1 bar 
R Hanode,Max 100 % — 
R Hanode, Min 0% = 
R Heathode, Max 100 % = 
R Heathode,Min 0% = 
Peathode,Max 3 bar 
Pcathode,Min 1 bar 
to 0 h 
£10,000 10,000 h 


365 


360 


355 


350 


345 


Temperature (K) 


340 — — Stack temperature at strategy 3 


— : Stack temperature at strategy 2 


335 — Stack temperature at strategy 1 








0 1000 2000 3000 4000 5000 «6000 7000 38000 9000 10000 
Hour 


Fig. 5 Temperature behavior in 10,000 h of fuel cell operation time in different strategies 


It can be infer from Fig. 6, that the anode and cathode optimum operating pressure 
are the same for all objective functions. It is also found that the cathode optimum 
operating pressure is higher due to water management considerations. 

Figure 7 presents the optimum humidity behavior in the three strategies. Inthe first 
objective functionsaturated situation is suggested. This is mainly due to maximum 
membrane ionic conductivity in saturated membrane. In the second objective func- 
tion, to minimize degradation due to hydrogen cross over lower relative humidity at 
anode and cathode is proposed. In other words, lower humidity causes less reaction 
on catalyst and causes lower catalyst degradation [15]. In third objective function the 
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Fig. 6 Pressure behavior in 10,000 h of fuel cell operation time in different strategies 


Relative humidity 





0 1000 32000 3000 4000 5000 6000 7000 $8000 9000 10000 
Time (hr) 


— Anode relative humidity at strategy 2 — -Cathode relative humidity at strategy 2 
—— Cathode relative humidity at strategy 1 = - -Anode relative humidity at strategy 3 





— = Cathode relative humidity at strategy 3 — — Anode relative humidity at strategy 1 


Fig. 7 Humidity behavior in 10,000 h of fuel cell operation time in different strategies 


optimum humidity is continuously decreased through time. It indicates that in long 
term applications the optimum humidity should be lower than saturated condition to 
avoid significant performance degradation. 

Figure 8 shows the behavior of fuel cell maximum power density through time. In 
the first objective function, power density is the most at initial time, but degraded dras- 
tically over time. In the second objective function regarding considered degradation, 
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Fig. 8 Power density in fuel cell operational life time in different strategies 


= 


Table 4 Life expectancy of three strategies at specific degradation percentage 


Membrane Life expectancy Catalyst layer Life expectancy 
degradation (%) degradation (%) 
Strategy | Strategy 2 Strategy 3 Strategy 1 Strategy 2 Strategy 3 
8 5,100 6,200 5,600 I5 100 1,500 200 
10 6,500 8,000 7,000 20 280 4,700 700 


12 8,100 10,000 — 8,900 25 1,400 10,000 2,200 


the maximum power density, is lower at initial time. However, the total energy over 
time is higher in comparison with the first objective function. In the third objective 
function, the maximum power density is evaluated in every desired time. Therefore, 
the total energy over time is much higher than the first and second objective functions. 
As the result, the total energy of system for third objective function is estimated to 
be 5.64 kW hr/cm?, which means 3.1696 higher in comparison with first objective 
function. 

It would be of interest if the optimization model is formulated to predict the life 
expectancy of fuel cell. To reach particular membrane and catalyst degradation, life 
expectancy of fuel cell at each strategy is presented in Table 4. 


6 Conclusion 


In this paper, a degradation model is developed to predict the effects of degradation on 
PEM fuel cell performance. The results indicate that through fuel cell lifetime, loss of 
catalyst active surface area is the main degradation mechanism, which directly affects 
voltage and power density loss. Moreover, operational conditions have a significant 
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impact on degradation. In the present work, three optimization strategies for maximiz- 
ing fuel cell electricity generation through lifetime by controlling operating conditions 
were compared. The objective function which optimizes the operating conditions in 
each time interval shows higher produced total energy. The optimal operating condi- 
tions are expressed as a function of time. It is found that by continuously optimizing the 
operating conditions, total energy generation of the system will increase up to 3.16 96. 
This study focused only on degradation mechanisms occurs atstationary applications 
of PEM fuel cell. In future works, a detailed comprehensive analysis of other degra- 
dation mechanisms will be discussed. 


Appendix: Nomenclature 


Model variables and constants 


CRef 
F 
kı 
k2 


t 

Cla 
2a 
ac 
Ue 

ul 
+max 
Au, 
DPt : PPtO 
ÓpiQ 
Opt 
OPtO 
NOH 
NH2 

F 

Q 

Ky 

K2 
PM, H2 
A PH? 


.H2 
lcross 


Fion 
tm 


A 


mol m^? 


C equiv! 
mol m~2s~—! 
mol m~25~—! 
kg 

kg mol! 

kg mol! 
equiv mol! 


equiv mol! 


J (mol K)! 
S 
Unitless 
Unitless 
Unitless 
V 

V 

mol m` 
J mol! 
Kg m^? 
Unitless 
J mol? 
J mol? 
mol s7! 
mol s7! 
Unitless 
J mol”! 
mol 
mol cm - 


2 


2 


Acmcm ?bar ! 


bar 

Mol cm ?57! 
Q 

cm 

Unitless 


Concentration refrence value (1 M = 1,000 mol m? ) 

Faraday's constant (96.485 C equiv-!) 

Pt dissolution reaction rate constant 

Pt oxidation reaction rate constant 

Total Pt mass in cathode 

Pt molecular weight (0.1951 kg mol ^ l) 

PtO molecular weight (0.2110 kg mol!) 

Number of equivalent electrons Transferred per mole of Pt 
dissoluted (2 equiv mol- 1) 

Number of equivalent electrons Transferred per mole of Pt 
oxidized (2 equiv mol~!) 

Universal gas constant (8.314 J mol-! K-!) 

Time 

Cathodic transfer coefficient of Pt dissolution reaction 

Anodic transfer coefficient of Pt oxidation reaction 

Cathodic transfer coefficient of Pt oxidation reaction 

Standard thermodynamic potential of Pt dissolution (1.188 V) 

Standard thermodynamic potential of Pt oxidation (0.980 V) 

Maximum oxide at Pt surface 2.18e—> (atom ratio of PtO = 1:1) 

Fitted chemical potential shift for Pt 

Pt metal density (21.090 kg m), PtO density (14.100 Kg m~?) 

fraction of platinum surface covered by PtO 

Surface tension of Pt nanoparticles (2.37 J m 2) 

Surface tension of PtO nanoparticles (1.00 J m 2) 

Rate at which radicals are produced 

Flux of hydrogen permeating through the membrane 

Fraction of remaining electrolyte 

PtO-PtO interaction coefficient 

Proportionality constant (1.8e—4 mol! ) 

Proportionality constant (1.8e—4 mol! ) 

Permeability of the membrane 

Partial pressure difference across the membrane 

Hə Crossover rate 

Membrane ionic resistance 

Membrane initial thickness (1 5e? cm) 

Membrane water content 
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a Unitless | Membrane relative humidity 

F cm Diffusion layer thickness (7.5e74 cm) 

Od cmQ-! Diffusion layer electronic conductivity 

I A cm ? Output current density 

in Acm ? Internal current density related to internal current losses 

i A cm~? Exchange current density related to activation losses 

P bar Fuel cell backpressure 

RH Unitless Relative humidity 

i], A cm~? Limiting current density 

T K Fuel cell temperature 

Py bar Saturated vapor pressure at humidifier temperature 

Psat bar Saturated pressure 

Vic V Fuel cell output voltage 

Vrev V Reversible voltage of fuel cell 

Volum V Ohmic loss 

V CrossOver V Activation and cross over loss 

VMassTransfer V Concentration loss 

p : bar H» Partial pressure 

piis bar Fuel cell anode inlet total pressure 

p Fii bar Hp Partial pressure in the inlet stream 

Po», Pg? bar Partial pressure of oxygen, partial pressure of hydrogen 
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